The cosmic 6 Li and 7 Li problems and BBN with long-lived charged massive particles 
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Charged massive particles (CHAMPs), when present during the Big Bang nucleosynthesis (BBN) 
era, may significantly alter the synthesis of light elements when compared to a standard BBN 
scenario. This is due to the formation of bound states with nuclei. This paper presents a detailed 
numerical and analytical analysis of such CHAMP BBN. All reactions important for predicting light- 
element yields are calculated within the Born approximation. Three priorly neglected effects are 
treated in detail: (a) photodestruction of bound states due to electromagnetic cascades induced by 
the CHAMP decay, (b) late-time efficient destruction/production of 2 H, 6 Li, and 7 Li due to reactions 
on charge Z = 1 nuclei bound to CHAMPs, and (c) CHAMP exchange between nuclei. Each of 
these effects may induce orders-of-magnitude changes in the final abundance yields. The study 
focusses on the impact of CHAMPs on a possible simultaneous solution of the 6 Li and 7 Li problems. 
It is shown that a priorly suggested simultaneous solution of the 6 Li and 7 Li problems for a relic 
decaying at t x ~ 1000 sec is only very weakly dependent on the relic being neutral or charged, unless 
its hadronic branching ratio is Bh <C 10 -4 very small. By use of a Monte-Carlo analysis it is shown 
that within CHAMP BBN the existence of further parameter space for a simultaneous solution of 
the 6 Li and 7 Li problem for long decay times t x ^ 10 6 sec seems possible but fairly unlikley. 



I. INTRODUCTION 

Big Bang nucleosynthesis (BBN) is one of the stan- 
dard pillars of modern cosmology. In its simplest version, 
reduced to a model with only one parameter, i.e. the 
contribution of baryons to the critical density, fl^h 2 w 
0.0224 standard BBN predicted and observationally 
inferred primordial light element abundances are very 
close. This holds particularly true for 2 H, and with some- 
what less confidence also for 4 He. However, when the 
A > 4 elements are considered agreement is less convinc- 
ing. The observationally inferred 7 Li/H ratio is about 
a factor three smaller than that predicted in SBBN Q. 
Moreover, 6 Li which is known to only be synthesized at 
the level 6 Li/H ~ 10~ 15 - 10~ 14 during SBBN has been 
recently observed in about a dozen metal-poor halo stars 
with abundance 6 Li/H - 3-5 x 10~ 12 0,[| . It is tantaliz- 
ing that these observations indicate a plateau-structure, 
similiar to that observed in 7 Li, i.e. 6 Li abundance inde- 
pendent of metallicity of the star, for stars at the lowest 
metallicities. A 6 Li plateau, should point to a pregalac- 
tic or primordial origin of this isotope, since the 6 Li had 
already been in place before stars produced metallicity 
(and cosmic rays). However, it is cautioned that fairly 
uncertain stellar pre-main-sequence (PMS) destruction 
of 6 Li could contrive to give an apparent plateau [5J. 

7 Li (as well as 6 Li) are observed in the atmospheres of 
metal-poor halo stars. When transported to the hotter 
interior of the star, by either convection or turbulence, 
both isotopes may be destroyed. It is thus possible that 
atmospheric 7 Li has been depleted by some factor though 
standard stellar models do not forsee this. A number of 
groups have recently re-studied this possibility [1, 0, HI • 
Postulating stellar turbulence with a parametrised mag- 
nitude, but of unknown origin, Korn et al Q claim that a 
star-to-star homogeneous factor 1.95 depletion is possible 
and even favorable when observations of the metal-poor 



globular cluster NGC6397 are considered. If true, the 
remaining factor ~ 1.5 could be either due to systematic 
errors in the effective stellar temperature calibration or 
due to an overestimate of the SBBN predicted 7 Li abun- 
dance due to systematic errors in nuclear reaction data. 
Concerning the second possibility, a recent remeasure- 
ment of the key 7 Li producing reaction ( 3 He(a, 7) 7 Be) 
seems to rather indicate a slight underestimate of the 
synthesized 7 Li |9j. 

6 Li is known to be produced by spallation (p + 
CNO — > LiBcB) and fusion (a + a — > Li) reactions by 
standard cosmic ray primaries scattering off nucleons and 
nuclei in the intergalactic medium [10( . Though this pro- 
cess may explain the observed 6 Li at solar metallicity, it is 
clear, however, that it falls short by a large factor (~ 50) 
to explain the 6 Li observed at low metallicity. Similiar 
holds true for putative cosmic ray populations due to 
shocks developed during structure formation [11]. In or- 
der to produce 6 Li/H~ 5 x 10~ 12 an earl y co smic ray pop- 
ulation of ~ lOOeV/nucleon is required Most candi- 
date sources fall short of this. The few viable remaining 
sources are due to accretion on the central Galactic black 
hole, albeit with an efficiency a factor 10 4 larger than 
that presently observed, or due to a significant fraction 
~ 0.1 of all baryons forming supermassive stars (and cos- 
mic rays) [12| . It may also be that our galaxy was host 
to a radio-loud quasar some time ago [la ]. The ener- 
getic problem becomes even exaggerated when likely 6 Li 
destruction during the stellar PMS [f| and putative 6 Li 
destruction during the stellar main-sequence phases are 
considered, possibly solving the 7 Li discrepancy. Finally, 
it has also been suggested that the 6 Li may result in situ 
from production by solar flares within the first billion 
of years of the star's life [13]. Though this seems pos- 
sible, it is hard to evaluate if a sufficient fraction of the 
freshly synthesized 6 Li falls back into the stellar atmo- 
sphere rather than being expelled by the solar wind. 
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It is entirely possible that the 7 Li and 6 Li anomalies 
are signs of physics beyond the standard model possi- 
bly connected to the quest for the cosmic dark matter. 
Even very small non-thermal perturbations in the early 
Universe may lead to a significant and observable 6 Li 
abundance, without overly perturbing other light ele- 
ments. It had thus been suggested that an anomalous 
high 6 Li abundance is due to non-thermal nuclear reac- 
tions (i.e. 3 H (a, n) 6 Li ) induced by the late-time 
t > 10 7 s electromagnetic [TEEH or hadronic [lTj decay of 
a relic particle, as for example the gravitino. 6 Li in abun- 
dance as observed in old stars may also be synthesized 
due to residual dark matter annihilations during the BBN 
epoch [TH. In particular, a standard thermal freeze-out 
process of weak scale particle dark matter (such as su- 
persymmetric neutralinos) is concommitant with the pro- 
duction of 6 Li in the right amount, given the dark matter 
mass falls in the range 20 < m x < 90 GeV, and annihi- 
lation is to a significant fraction hadronic and s-wave. 
Concerning a solution to the 7 Li problem, early attempts 
utilising the electromagnetic decay of a relic and the in- 
duced 7 Be photodisintegration [l9| ( 7 Li is mostly syn- 
thesized as 7 Be, which later on electron-captures) have 
not proven viable due to una ccep table perturbations in 
the 2 H/H and 3 He/ 2 H ratios [20]. However, it has been 
shown that the hadronic decay of a relic during BBN, 
and the induced excessive neutron abundance may pre- 
maturely convert 7 Be to 7 Li which is then destroyed by 
proton capture. When fl x Bh ~ 1 — 5 x 10~ 4 , where Bh is 
the hadronic branching ratio, a factor 2 — 4 destruction 
of 7 Li results [2l[. For relic decay times w 1000 s, it is 
moreover possible to synthesize all the observed 6 Li by 
non-thermal nuclear fusion. This has been the first, and 
so far only, known simultaneous solution to the 6 Li and 
7 Li problems. It is noted that such a decay also leads to 
a possibly problematic 30% - 50% increase in the synthe- 
sized 2 H/ 1 H ratio. 

Within the context of minimal supersymmetric exten- 
sions of the standard model of particle physics, a simul- 
taneous solution is nicely realised, either by heavy grav- 
itino decay, or in the case that gravitinos are the lightest 
supersymmetric particles (LSPs) by the supersymmetric 
partner of the tau-lepton (the stau) decaying into grav- 
itinos [2l|. In the second scenario, an added benefit is 
that for the right parameters to solve the 6 Li and 7 Li 
problems, TeV staus left over from a thermal freeze-out 
at higher temperature, and decaying at t x k 1000 s into 
50 — 100 GeV gravitinos produce naturally about the right 
amount of gravitinos to explain the dark matter and of 
a warmness interesting to the formation of large scale 
structure formation [22j |. Unfortunately, staus of mass 
1 TeV are too heavy to be discovered at the LHC. 

Recently, it has been realised that the existence of elec- 
trically charged massive particles (CHAMPs) during the 
BBN epoch may lead to modifications of the synthesis 
of light elements [H, [24], HE] beyond those simply due to 
their decay. Since for gravitino LSPs, the next-to-LSP 
(NLSP) is long-lived and in about half of the supersym- 



metric parameter space it is the electrically charged stau, 
such effects arc important to consider. Other metastable 
charged relic particl es p ossibly existing during BBN have 
been also proposed [26| . Modifications to BBN occur due 
to the formation of electrically bound states between the 
negatively charged CHAMPs and the positively charged 
nuclei. The realization that (meta)-stable weak-scale 
mass charged particles enter into bound states during 
and after BBN had already been made in the late eight- 
ies [13, [H, [29[ , when the possibility of charged dark mat- 
ter was analyzed. Nevertheless, the influence of bound 
states on BBN had not been much discussed. 

In this paper results of the up-to-now most detailed 
calculations of BBN nucleosynthesis in the presence of 
decaying negatively charged particles are presented. The 
analysis attempts to reveal all key processes important 
for a reliable prediction of light element yields, thereby 
revealing, heretofore neglected effects, which make orders 
of magnitude changes in the predicted BBN yields for 
much of the parameter space. These changes are found 
mostly for late decaying t x > 10 6 s CHAMPs. The aim 
of the paper is to analyze the potential of bound-state 
nucleosynthesis to solve the cosmic 7 Li and 6 Li problems. 

The outline of the paper is as follows. In Section 2 
a discussion/analysis of all prior ly suggested solutions 
to the 7 Li problem within bound-state nucleosynthesis 
is presented, whereas in Section 3 details of the present 
calculations are given. In Section 4 it is shown that BBN 
continues to very low temperatures T -C 1 keV in the 
presence of bound states. Section 5 shows that bound 
states are efficiently photodisintegrated already at high 
temperature due to the decay of the relic. Section 6 
stresses the importance of CHAMP transfer reactions at 
late times. Finally In Section 7 possible further solutions 
to the 6 Li and 7 Li problems for late-decaying CHAMPs 
t x ^ 10 6 sec are discussed. Section 8 draws the conlusions. 
An appendix gives some detail on the determination of 
reaction rates in the Born approximation. 



II. BOUND-STATE BBN AND PRIOR 
SUGGESTED SOLUTIONS TO THE 7 LI 
PROBLEM 

Modifications to BBN occur due to the formation of 
electrically bound states between the negatively charged 
CHAMPs and the positively charged nuclei. Since bound 
state binding energies may be appreciable (cf. Table 
1), a significant fraction of Be may be captured by 
CHAMPs at temperatures as high as T < 30 keV, whereas 
the same occurs at T < 10 keV for 4 He. This may be seen 
in Fig.[TJ which shows the fractions j\ — n^.x-) K\ of 
7 Be, 7 Li, 6 Li, and 4 He locked up within bound states. On 
first sight, the most important effect of bound states dur- 
ing BBN is a reduction of the Coulomb barrier [23|, [24j . 
Nevertheless, since SBBN is essentially finished at T m 
10 keV, Coulomb barrier modifications of reactions rates 
involving He should be hardly important (even though, 
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ad hoc, speculated otherwise in Ref. [24]]). However, 
as shown by Pospelov 23] there is a non-trivial cat- 
alytic effect on reactions involving photons in the final 
state. SBBN reaction rates involving dipole radiation 
(El; e.g. 3 He( 4 He, 7) 7 Be) scale as A~ 3 , whereas reac- 
tion rates forbidden at the dipole approximation but al- 
lowed at quadrupole (E2; e.g. 2 H( 4 He, 7) 6 Li) scale as 
A~ 5 , where A 7 is the wavelength of the emitted pho- 
ton. This, in both cases is around ~ 130 fm. In the 
presence of a 4 He-CHAMP bound state the reaction may 
proceed photonless (e.g., 2 H( 4 Hc— X~, A~) 6 Li) and A 7 
is approximately replaced by the Bohr radius a4 Ho of the 
4 He-CHAMP bound system. Since a4 Hc « 4.8 fm (cf. 
Table 1) very large enhancement factors of 7 x 10 7 and 
3 x 10 5 [13, IS^ for the S-factors of the 2 H + 4 He, and 
3 He + 4 He reactions, respectively, have been estimated. 
A recent more detailed three-body nuclear reaction calcu- 
lation of the 2 H + 4 He reaction, has reduced this estimate 
by a factor ~ 10 3l|. Such large enhancement factors are 
important as they lead to excessive 6 Li (and 7 Li) pro- 
duction for any weak scale charged particles which are 
sufficiently long-lived t x > 4 x 10 3 s, unless fix *3x 10~ 6 . 
They have thus been utilised to place a stringent upper 
limit on the reheat temperature in the early Universe 
T< 10 7 GeV in the case when the supersymmetric grav- 
itino exists and when it is the LSP [321] . Nevertheless, it 
seems somewhat premature to set such upper limits, as 
the BBN with charged long-lived particles for decay times 
tx ^ 10 6 s had priorly not been investigated (cf. Section 



The putative existence of bound states during BBN 
has also led to a flood of claims of possible solutions to 
the 7 Li and/or 6 Li anomalies. In Ref. [24[ it was real- 
ized that significant fractions of the 7 Be and 7 Li isotopes 
are within bound states during BBN. This has lead the 
authors to arbritrarily enhance certain reactions rates 
involving mass-7 element destruction processes by large 
factors, leading to the claim that the existence of bound 
states may solve the Li overproduction problem. How- 
ever, these claims are, up to now, unfounded [33[ (see 
also below) . In Ref. [25] it was noted that during the de- 
cay of A - , when residing in a bound state with He, the 
4 He nucleus could break up. The resultant energetic 3 H 
and 3 He could then fuse on 4 He to produce 6 Li, in a sim- 
ilar to what had been proposed in [TBI, [l7| ■ Though the 
suggestion is correct, the authors calculate the break-up 
probability to be very small (cf. also Ref. [35] ) , such that 
the 6 Li synthesis by catalytic 2 H ( 4 He-A~, A~) 6 Li is by 
far dominant. The analysis of Ref. 30 (and Ref. (34|) 
essentially confirms the simultaneous solutions to the 6 Li 
and 7 Li problems as given in Ref. [2ll . |22| | , e ven when 
bound state effects are included. In Ref. [3(| the case 
of almost degenerate NLSP staus f and LSP neutrali- 
nos x na s been considered. Here mass splittings smaller 
than 5m = rrif — < 1 GeV have been assumed. In 
this region of f~x parameter space, motivated by the 
well-known f-\ coannihilation region for neutralino dark 
matter, the stau is relatively long-lived due to final phase 
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FIG. 1: Bound state fractions j\ = n^ N . x -)/nN* of nuclei Ni 
bound to CHAMP X~ as a function of temperature T, for a 
model with Mx = 100 GeV and Qxh 2 =0.1 (corresponding 
to a CHAMP-to-baryon ratio Y x - = 4.26 x 10" 2 /2). Shown 
are fi for 7 Be solid (red), 7 Li long-dashed (green), 6 Li short- 
dashed (blue), and 4 He dotted (purple), respectively. Nuclear 
destruction of bound states results in a behaviour of ff dif- 
ferent than that expected from simple estimates by the Saha 
equation. This is particularly seen in /f for 7 Li due to the 
iHfLi-X", A") 4 He + 4 He reaction. 



space supression of the decay. It is claimed, that the 7 Li 
overproduction problem may be solved by internal con- 
version of staus in bound states with 7 Be, to neutralinos, 
e.g. (f- 7 Bc)^ x + ^r+ 7 Li and the subsequent destruc- 
tion of 7 Li by protons. It is argued that solutions to the 
7 Li problem may be found for 8m < 100 MeV even for the 
smallest abundances of staus. A more detailed analysis 
of the 7 Be-bound state fraction via the Boltzmann equa- 
tion shows, however, that only a very small fraction of 
7 Be are within bound states, thus making modifications 
of the 7 Li abundance at low stau-density negligible. At 
larger stau-densities some effect may result. 

Ref. 35] make the interesting suggestion that the 7 Li 
problem could be solved by catalytic conversion of 7 Be 
via ( 7 Be-X~)(p, 7)( 8 B-A~) and the subsequent beta- 
decay of the 8 B -^ 8 Be+e + + v e nucleus. This reac- 
tion would mostly occur via a p- 7 Bc resonance in the 
8 B nucleus which, in the absence of bound states lies at 
769.5keV relative to the p- 7 Be continuum. The catalysm 
in the reaction would then occur by a shifting of the 
resonance to « 167 keV relative to the p-( 7 Be-A _ ) con- 
tinuum since the ( 8 B— A~) bound state binding en- 
ergy (Ea BX - ~ 2.0 MeV) is larger than that of 7 Be 
{Etq b x- ~ 1-39 MeV), making the resonance available 
at only slightly supra-thermal energies. Moreover, apart 
from the decrease in the resonance energy they also 
deduce a factor ~ 10 3 larger reaction rate coefficient. 
Adopting their calculated rates for 7 Be-A~ bound state 
formation and the ( 7 Be-X~)(p, 7)( 8 B-A~) reaction, I 
partially confirm this effect by full numerical analysis. 
For example, for tx — 1-5 x 10 3 s and the number ratio 
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of X s to baryons Yx ~ 0.2, I find a reduction of the 
7 Li abundance by 33%, and a 6 Li/H ratio of 2 x 10~ u . 
However, the effect is not as strong as initially imag- 
ined, since by the reciprocity theorem the inverse rate 
is also enhanced. The inverse rate l/ri nv is thus around 
10 3 times larger at T ss 32.5 keV than the beta decay 
rate I/173 of 8 B (half-time of 770 ms), converting 8 B-X~ 
rapidly back to 7 Be-X~ + p, before 8 B can beta-decay. 
The effect is therefore essentially absent at early times 
(i.e. small tx)- Nevertheless, the inverse rate quickly 
drops below the beta decay rate (i.e. Tp/T urv « 0.1 at 
T as 24.1 keV). For the same parameters as above, I still 
find a 14% reduction of the final 7 Li. This drops to 7% , 
2% for tx = 10 3 and 7 x 10 2 s, respectively. 

It is interesting to know if the solution of the lithium 
problems proposed in Ref. [2l| is changed when the de- 
caying relic is charged, such as the stau. In Fig. [J] the 
parameter space solving either the 7 Li problem, or both 
the 6 Li and Li problems, is shown. The upper panel 
shows results for a charged relic and the lower panel for 
a neutral relic. Here observational limits as discussed in 
Ref. (3?| have been applied and the 6 Li, 7 Li problems 
are assumed to be reconciled with observational data for 
6 Li/ 7 Li >0.03 and la/^H < 2.5 x lO" 10 . The assumed 
parameters of the model are a hadronic branching ratio 
Bh = 10~ 4 and relic mass Mx — 1 TeV. It is seen that 
even at Bh as small as 10 -4 the 7 Li-solving region is es- 
sentially unmodified, whereas some changes are observed 
in the 6 Li and 7 Li solving regions. These latter are mostly 
due to excessive 6 Li production when the relic is charged, 
disallowing some of the larger life times tx ^ 2 x 10 4 s. 
Bound state effects are nevertheless important when the 
hadronic branching ratio is very small. This may be seen 
in the lowest panel of Fig. [21 where Bh = has been as- 
sumed. When only bound state effects are operative, the 
2 H/ 1 H-ratio is essentially unmodified. This is in contrast 
to the solution of the lithium problems with a hadronic 
decay, as seen by the dotted (blue) lines in the upper 
two panels, beyond which 2 H/ 1 H is larger than 4 x 10~ 5 . 
It is intriguing that both processes, hadronic decay and 
bound state effects, have the same preferred tx for a 
simultaneous solution of the lithium problems. 



III. DETAILED BOUND-STATE BBN 
CALCULATIONS 

The calculations presented here attempt to take proper 
account of the influence of singly bound states on the nu- 
cleosynthesis for elements with nucleon number A < 7. 
Heavier elements as well as the formation of molecules, 
such as (X~ — 4 He — X^) 1 are not considered. All effects 
of electromagnetic and hadronic cascade nucleosynthe- 
sis are included and treated as presented in Ref. [37| . 
The fractions of inividual nuclei i in bound states ff = 
n (N i x-)/ r! ' t N t are computed by full numerical integration 
of the Boltzman equation. This is required since esti- 
mates by the Saha equation are only very approximative, 



TABLE I: Nucleus, energy of bound state, approximative 
Bohr radius of bound state as [38|, and adopted root-mean- 
square charge radius for nucleus 



nucleus 


rjb (Ke\ ) 


~ as (fm) 


/„2\ 1/2/ , \ 

(r ) c (fm) 


J H 


24.97 


28.8 


0.895 


2 H 


49.5 


14.4 


1.3 


3 H 


72.6 


9.6 


1.7 


3 He 


269 


5.2 


1.951 


4 He 


349.6 


4.8 


1.673 


6 Li 


842.5 


2.1 


2.37 


7 Li 


897.6 


1.9 


2.50 


7 Be 


1385 


1.5 


2.50 



due to the relatively early freeze-out of the CHAMP- 
nuclei recombination process [24[. Except of the recom- 
bination rate of on 7 Be, which is taken from Ref. [35j . 
all other recombination rates are computed by a numer- 
ical integration of the Schroedinger equation. This may 
make difference up to a factor two in f\ since the recom- 
bination rates as given in Ref. [24j only apply asymptot- 
ically at low temperature T . Bound state wave functions 
and bound-state energies are also computed by an inte- 
gration of the Schroedinger equation, assuming realistic 
charge radii for the nucleus as measured by experiment. 
The reader is referred to Table 1, for some of the bound 
state properties. Finally, it is important, to take into ac- 
count the nuclear destruction of bound states. Nuclear 
rates are very fast at early times, and for reaction which 
are sufficiently exothermic, the electric bound between 
the final nucleus (nuclei) ought to be destroyed [39] . This 
often changes j\ by orders of magnitude. 

A proper evaluation of BBN yields with bound states 
is only possible when somewhat realistic nuclear reaction 
rates for nuclei within bound states are present. With the 
exception of the reaction 2 H( 4 He— X~, X~) 6 Li a more 
detailed evaluation of such reactions had been absent of 
the literature so far. Improving over simple scaling re- 
lations [23], [3(j seems important also, since nuclear reac- 
tions including bound states contain three quantities of 
similar magnitude, 0,3 the Bohr ratius, a nuc i the nuclear 
nuclear radius, and kf the momentum of the outgoing 
nucleus. All three quantities are in the several Fermi 
range, thus leading potentially to important cancellation 
effects. More importantly, estimates via simple scaling 
relations adopt the Born approximation, which is known 
to fail at low energies and strong perturbations [40] . This 
is essentially the case for all reactions of importance to 
bound-state BBN. The failure of the Born approxima- 
tion had been seen, for example, by the reduction of the 
2 H( 4 He— X~ , X~) 6 Li rate by a factor ~ 10, when a more 
detailed evaluation [3l[ is compared to a simple scaling 
result. 

I have identified all key reactions in bound-state BBN. 
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FIG. 2: Parameter space in the relic particle-to-baryon ra- 
tio Yx and relic particle decay time tx which may resolve 
either the 7 Li problem (green - light) or both, the 7 Li and 
6 Li problems (red - dark). The panels show, from top to bot- 
tom: (a) a charged relic with Bh = 10 -4 , (b) a neutral relic 
with Bh = 10~ 4 , and (c) a charged relic with Bh = 0. All 
three panels assume a mass M x = 1 TeV for the relic. By 
comparison of the green (lighter) areas it is seen that bound- 
state effects on 7 Li, as suggested in Ref. [35), do not have a 
very big impact for relic hadronic branching ratios Bh 1CP 4 . 
The adopted abundance limits are: 2 H/ 1 H< 5.3 x 1CP 5 , 
7 Li/ 1 H< 2.5 x 1(T 10 , 6 Li/ 7 Li< 0.66, and 6 Li/ 7 Li> 0.03 to 
solve the 6 Li problem. Above the dotted lines the 2 H/ 1 H 
ratio exceeds a value of 4 x 10 -5 . 



These are shown in Table 2. It is completely beyond the 
scope of the present paper to evaluate all these reaction 
rates more properly, i.e. beyond the Born approxima- 
tion, a task which is formidable in particular when the 
important CHAMP-exchange reactions (cf. Section 6) 
are also considered. For the 2 H( 4 He— X~, AT~) 6 Li pro- 
cess the rate as given by Ref. 31] was adopted. For other 
reactions, as a starting point, I have thus nevertheless, 
evaluated rates in the Born approximation. These rates 
will serve as benchmarks later on. For details concern- 
ing these calculations the reader is referred to Appendix 
A. Results for the in such a way obtained S-factors are 
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FIG. 3: Nuclear reaction S'(£')-factors as function of 
energy computed in the present analysis. The most 
important 5*-factors for nuclear reactions involving the 
( 4 He-X") bound state are shown: 2 H( 4 He-X" , A"-) 6 Li 
solid (red), 3 H( 4 He-X",X-) 7 Li short-dashed (blue), and 
3 He( 4 He-X~,X~) 7 Be dotted (purple), respectively. The 
dashed-dotted (light-blue) line shows the result of a recent 
evaluation [3l[ of 2 H( 4 He— X - , X~) 6 Li, whereas the long- 
dashed (green) line shows the result for the same reaction 
computed in this paper when the / = 1 and I — 2 contribu- 
tions are neglected. 



TABLE II: Assumed properties for the calculation of nuclear 
reactions with one nuclei in a bound state. The columns show: 
Reaction, S-factor for the SBBN reaction in MeV barn, angu- 
lar momentum for the (AB) = C final bound nucleus, and the 
multipoles for the initial Coulomb wave which are included in 
the calculation. 



No. 


(AX) + B ->• C + X 




Ic 


l Coul 


1 


( 4 He-A~) + 2 H^ 6 Li + X- 


10" 8 





0,1,2 


2 


( 4 Re-X~) + 3 H^ 7 Li + X~ 


8 x 10~ 5 


1 


0,1 


3 


( 4 He-A~) + 3 He^ 7 Be + X~ 


4 x 10 -4 


1 


0,1 


4 


CR-X-) + 6 Li^ 7 Bc + X~ 


10~ 4 


1 


0,1 


5 


CR-X~) + s Lh 4 He + 3 He + X~ 


3 






6 


( 1 H-X~) + 7 Li^ ( 8 Be-X~) + 7 


10 -3 


1 


0,1 


7 


CR-X-) + 7 Be^ 8 B + X~ 


3 x 10~ 5 


1 


0,1 


8 


( 2 H-A^) + 4 He^ 6 Li + X~ 


10" 8 





0,1,2 


9 


( 3 H-A") + 4 He^ 7 Li + X- 


8 x 10 -5 


1 


0,1 



shown in Fig. [3] and Fig.|H respectively. 



IV. LATE-TIME BOUND-STATE BIG BANG 
NUCLEOSYNTHESIS 

The reader may have noted that Table 2 also includes 
reactions with bound states on elements 1 H, 2 H, 3 H with 
only one charge number Z — 1. In fact, such reactions are 
extremely important at low temperatures T < 3, 2, 1 keV 
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FIG. 4: Nuclear reaction S , (i5)-factors as function of en- 
ergy computed in the present analysis. The most impor- 
tant 5-factors for nuclear reactions involving bound states 
with Z = 1 nuclei are shown: (from top to bottom 
at the highest energies) 6 Li( 1 H-X _ , X^) 4 He+ 3 He double- 
dotted (black), 7 Be( L 'H-X~,X~) 8 B long-dashed (green), 
4 He( 3 H-X~,X~) 7 Li dotted (purple), 4 He( 2 H-X-,X-) 6 Li 
dash-dotted (light-blue), 7 Li( 1 H-X~, 7 )( 8 Be - X") short- 
dashed (blue), and 6 Li( 1 H-X~, X") 7 Be solid (red). 



when one after the other, non-negligible fraction of 3 H, 
2 H, and 1 H enter into bound states. This may be seen 
in Fig. [5j It is noted here, that a possible impact of 
such reactions has been pointed out before [24[ , albeit in 
a very approximative way. It was not clear, a priori, if 
the Coulomb barrier between, for example, 1 H and 6 Li is 
sufficiently supressed in order to make reactions such as 
6 Li( 1 H-X", X~) 4 Ke + 3 He efficient enough to substan- 
tially reduce any priorly synthesized 6 Li. This is because, 
on first sight, Coulomb shielding of the proton could 
only be partial, due to the fairly extended Bohr radius 
as ~ 29 fm of the 1 H-X~ system. In Fig.[5]a I — spher- 
ical wave without any Coulomb repulsion, i.e. V c = 0, is 
compared to the spherical Coulomb wave functions be- 
tween the 6 Li and the 1 H-X~ bound state with I = and 
I = 1 initial angular momentum, respectively. It is seen 
that essentially no Coulomb supression exists. Rather, 
the incoming wave function of the 6 Li nuclei is even 
strongly enhanced at the center. This is not surprising, 
as by assumption, the X~ resides at the center, and due 
to the significant spread in the wave function of the pro- 
ton (clb ~ 29 fm) the effective proton charge density at 
the center is low. The Coulomb potential for the 6 Li nu- 
cleus is (/>6Li = — 3e 2 exp(— 2r/as)(l/r + 1/as), thus very 
attractive at the center and approaching zero at large 
distances. Nuclear reactions between such bound states 
and bare nuclei, are therefore not Coulomb supressed. It 
is rather conceivable, that Coulomb focussing occurs at 
low energies, even enhancing the reaction rates over the 
Vc = case. This may be observed in the S-factor for the 
(t-H-X-) + 6 Li^ 4 He + 3 He + X~ reaction as shown 
in Fig. |U It is noted here, that due to an anomously 
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FIG. 5: Bound state fractions f\ for 3 H (solid - red), 2 H 
(dashed - green), and 1 H (blue - dotted) as a function of tem- 
perature T. Adopted model parameters are as in Fig. [T] with 
a X decay time rx = 10 s. For illustrative purposes photo- 
disintegration of bound states due to X-decay (cf. Section 5) 
and X-exchange reactions (cf. Section 6) have not been taken 
into account. 



low 7 Li( 1 H— X~, X~) a Be rate found in the Born approx- 
imation the rate for 7 Li( 1 H— X~, 7)( 8 Be— X~) has been 
comuputed and utilised in the calculations. 

Thus, Z = 1 bound states at T w 1 keV behave 
almost as neutrons (with the exception that they are 
stable). Already very small fractions of these bound 
states induce therefore a second round of late-time nu- 
cleosynthesis, capable of destroying all the synthesized 
6 Li, 9 Be, and some of the 7 Li. This may be seen in Fig. [7] 
where the 6 Li/H, 7 Li/H, 7 Be/H, and 2 H/H ratios are 
shown for a CHAMP with fl x h 2 = 0.01, m x = 100 
GeV, and decay time tx — 10 10 s, where h is the dimen- 
sionless present-day Hubble parameter, and fix HU the 
fractional contribution of CHAMPs to the present crit- 
ical density would they not have decayed. Note, that 
this is easily converted to the CHAMP-to-baryon ratio 
Y x = (tt x h 2 /n b h 2 )(m p /m x ) which is Y x ~ 4.26 x 10~ 3 
for the adopted parameters. The calculations presented 
in Fig. [7] (as well as Figs. Q] and [5]) are performed un- 
der the assumption that the X decay is not associ- 
ated with any electromagnetic- or hadronic- energy re- 
lease and in the absence of X-exchange reactions (cf. 
Section 6). This is done to isolate the effects of the 
bound states. At early times, towards the end of con- 
ventional BBN, when a significant fraction of 4 He en- 
ters bound states, the reactions 2 H( 4 He— X~, A"~) 6 Li, 
3 H( 4 He-X-,AT-) 7 Li, and 3 He( 4 He-AT~, AT~) 7 Be, syn- 
thesize significant, and observationally completely unac- 
ceptable abundances of the A > 4 isotopes. However, 
when bound states of the Z = 1 elements form at T w 
1 keV, essentially all the synthesized 6 Li and 7 Be may be 
rapidly destroyed by the reactions 6 Li( 4 H— X~, AT~) 4 He 
+ 3 He and 7 Be( 1 H-X", X~) 8 B. The situation ap- 
pears different for 7 Li, due to the small estimate for 
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FIG. 6: Spherical Coulomb wave functions of a Li nuclei 
with energy E = IkeV in the electric field of the 1 H-X~ 
bound state, for s-wave (angular momentum 1 = 0- solid - 
red) and p-wave (angular momentum 1 = 1- dotted - blue) 
where X~ is at radius r = 0. For comparison the spherical 
wave function without any Coulomb barrier, i.e. Vc = 0, for 
s-wave, is also shown (dashed - green). It is seen that no 
significant Coulomb barrier supression of the wave function 
near the origin exists. Rather, both Coulomb wave functions 
are significantly enhanced at the center, due to the presence 
of X~ at r — 0. The oscillatory behaviour may lead to im- 
portant interference effects. Both, the 1 = and I = 1 initial 
states have significant contributions to the cross section. 



the 7 Li( 1 H-X-,X-) 8 Be and 7 Li( 1 H-X~, 7 ) 8 Be-X~ 
cross sections, implying that almost all initially syn- 
thesized 7 Li is left intact [42]. The abundance of 
7 Li/H is found at an observationally friendly 2.7 x 
10~ 10 . It is noted that 2 H is also destroyed, though 
to a much smaller degree, mostly by the reactions 
3 H( 2 H-X-,n) 4 He+X-, 3 He( 2 H-A-,p) 4 He+X-, and 
2 H( 3 H— X~, n) 4 He+A"~, and to a lesser degree by 
2 H( 1 H-X-,X-) 3 He. The reader is referred to Table 
3 concerning assumptions about the rate of these, and 
some other reactions involving only A < 4 elements. 
Furthermore, when regarding Fig. [7] in more detail, one 
also notes late-time production of 6 Li and 7 Be at some 
level due to the 4 He( 2 H-A~, X~) 6 U as well as the 
"Li^H-X-, X~) 7 Be reactions. 

It is thus premature to conclude, that extreme 6 Li 
overproduction, rules out the existence of CHAMPs with 
long life times [23|, |43| . The model shown above, at 
CHAMP densities many (five !) orders above those al- 
ready claimed to be ruled out by 6 Li overproduction is ob- 
servationally viable in all abundances. Constraints on the 
existence of CHAMPs in the early Universe could there- 
fore, in principle, be much milder for long X~ life times 
than initially predicted. Nevertheless, they is further im- 
portant physics entering the calculations discussed in the 
next two sections. 



TABLE III: Assumed enhancement factor of a number of nu- 
clear reactions between A < 4 nuclei involving bound states 
of Z = 1 nuclei. The Coulomb supression factor is assumed 
to be completely absent in these reactions. 



No. 


(AX) + B->C + X 


enhancement 


10 


2 H( 1 H-X-,X") 3 He 


1.25 x 10 2 




1 H( 2 H-X-,A") 3 He 




11 


3 H( 1 H-X-,A-) 4 He 


10.7 




1 H( 3 H-X-,A~) 4 He 




12 


2 H( 3 H-X-,n) 4 He+A- 


1 




3 H( 2 H-A"-,n) 4 Ffe+X- 




13 


3 He( 2 H-X-,p) 4 He+X- 


1 
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(solid - red), 7 Li/ 1 H (long-dashed 
dashed - blue), and 2 H/ 1 H (dotted 



Be/ 1 !! 
8 Li/ 1 H (short- 
for a CHAMP 

model with Mx = 100 GeV, Q. x h 2 = 0.01, and r x = 10 10 s. 
It is seen that large amounts of 6 Li and 7 Be synthesized at 
T rs 10 keV will be again destroyed at T w 1 keV. Neither 
effects due to electromagnetic and hadronic energy release 
during CHAMP decay nor charge exchange effects have been 
taken into account. 



V. PHOTODISINTEGRATION OF BOUND 
STATES BY THE DECAY OF THE CHAMPS 

There is another effect, heretofore overlooked, which 
may significantly reduce the in catalytic BBN at T w 
10 keV synthesized 6 Li (and 7 Li) abundance. CHAMP 
decays are typically accompanied by the injection of elec- 
tromagnetically interacting partcles, with total energy 
comprising often a large fraction of the X rest mass. 
It is well-known, that such particles (e~,e + , and 7's) 
induce a rapid cascade on the cosmic blackbody pho- 
tons, due to 77bb pair creation and inverse Compton 
scattering e + 7 processes, until the energy of any re- 
maining 7's is too low to further pair-produce, i.e. for 
' th ~ m 2 /22T w l^MeVfjyiOkeV)- 1 . It is 



E y < E, 
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seen, that this energy is above the binding energy of 
4 Hc-AT~ (and 1 H-X _ ) even at temperatures as high as 
T w 30keV, making possible the 4 He-X" and ( X H-X-) 
bound state photodisintcgration before any significant 
6 Li synthesis (destruction) has occured. In Fig [8] the 
resultant photon spectrum due to the injection of en- 
ergetic electromagnetically interacting particles at cos- 
mic epochs with temperature T — 10, 1, and 0.1 keV is 
shown. The shown spectrum E 7 dn 7 / dlrtEy is generated 
by a Monte-Carlo simulation taking account, not only of 
pair production and inverse Compton scattering, but 
also 77 scattering (important at high Ey<E t h, Bcthc- 
Heitler pair production 7 + p, 4 He — > p, 4 He + e~ + e+, 
Compton scattering of the produced e ± , as well as the 
important Thomson (Klein-Nishina) scattering of 7's on 
thermal electrons. It is based on the calculations pre- 
sented in Ref . [37[ , with the Thomson scattering process 
extended to energies as low as E 1 « 25keV, to account 
for 1 H-X~ destruction. 

Following secondary and tertiary, etc. generations of 
scattered photons to obtain the correct photon spectrum 
for the bound state destructions process is mandatory. 
For example, the injection of 1 TeV of electromagneti- 
cally interacting energy at T = 1 keV is associated with 
injection of N 7 w 3.3 x 10 6 primary photons with en- 
ergy E 1 > 25 keV, resulting from the initial cascade on 
the blackbody. When further interactions of these 7's 



are considered the number rises to AT- 



1.1 x 10 s . In 



other words, an injected photon takes about 30 inter- 
actions before dropping below the threshold for 4 H-AT~ 
photodisintegration. This exemplifies the importance of 
subsequent 7 interactions. In Fig.[8]one may note a "pile- 
up" of photons at low E 1 . This is due to the typical frac- 
tional loss of 7's in the Thomson regime E~ < m e being 
small, such that it takes several Thomson scatterings for 
a photon to have dropped below £ 7 < _E? H w 25 keV. A 
similiar pile-up does not exist at E 1 < S| He w 350 keV 
since during scatterings of 7's with energy E 7 ~ m e on 
electrons the 7's may loose a significant fraction of their 
energy. We thus expect the effects of photodisintegra- 
tion of bound states have a larger impact on the 4 H-AT~ 
bound state fraction than on that of 4 He-V~. This effect 
is not only due to the above, but also due to the photodis- 
integration cross section of 4 H-A"~, c7 H _x- being larger 
than the one for 4 He-V _ . Note that all calculations be- 
low, include numerically evaluated cross sections for the 
photodisintegration of all A < 7 nuclei bound states. 

In Fig. [5] the bound state fractions in two scenar- 
ios: (a) of 4 He for a model with Vtxh 2 — 0.1 and 
tx = 3 x 10 4 s (and electromagnetic decay), and (b) of 
4 H for VL x h 2 = 5x 10~ 3 and t x = 3 x 10 6 s, are shown 
in the same graph. Here the solid lines show /| H (/i b H ) 
when non-thermal bound state photodestruction is in- 
cluded, whereas the dotted lines show results when it is 
neglected. It is seen that realistic bound state fractions 
are significantly lower. In scenario (a) a 6 Li/H ratio ~ 10 
times lower results, compared to when photodestruction 
is neglected, whereas in scenario (b) the 6 Li/H ratio is 




E (MeV) 

FIG. 8: Resultant photon spectrum _E 7 dn 7 /d ln_E 7 due to 
electromagnetic energy injection at cosmic epochs with tem- 
peratures at T = lOkeV (solid -red), IkeV (dashed - green), 
and 0.1 keV (dotted - blue), respectively. The normalisa- 
tion of the spectrum is arbritrary. The fraction of photons 
with energy above the 4 He-X~ photodisintegration threshold 
E\ « 350 keV is w 1.9%, 3.7%, and 4.7% for temperatures 
T = 10, 1, and 0.1 keV, respectively. 



~ 100 times higher. Here case (b) is affected by a reduced 
efficiency of 6 Li( 1 H— X~, V~) 4 Hc + 3 He, whereas in case 
(a) the reaction 2 H( 4 He— X~, AT~) 6 Li is rendered less 
dominant. For sufficently high Qx, and when thermal 
photodisintegration is unimportant, the resultant bound 
state fraction may be estimated by a steady state between 
the recombination rate, i.e. {av) rec n<tYL C nx- and the 
photodisintegration rate, i.e. (crc) p hn(4 He _ x -)n 7 . Here 
rt4 Hc , n(4 Ho _ x -), n x -, and n 7 are free 4 He, bound 4 He, 
X~, and nonthermal photon number densities, respec- 
tively. The nonthermal photon number density n 7 may 
be obtained from n 7 « dnx /dt tthN^ where dnx/dt w 
nx /i~x before substantial decay, tth is the life time of 
photons against Thomson scattering (i.e. the typical sur- 
vival time), and N^ b is the typical number of photons per 
particle decay with energy above the photodisintegration 
threshold Eb (including secondary generations). This, 
for example at T = 1 keV, is approximately 4 x 10 6 and 
1 x 10 8 for 4 He and X H bound state photodisintegration, 
respectively, per 1 TeV of electromagnetically interacting 
energy injected into the plasma. It is thus found 



n 



rc( 4 He-X-) (cv)rcc T X 1 



He 



714 Hc ((Jc)ph T T h A^S b 



(1) 



It may be noted that this expression, which is valid only 
for large Y x > 10~ 2 is independent of the CHAMP-to- 
baryon ratio, but dependent on the CHAMP life time. 



VI. CHAMP EXCHANGE REACTIONS 



It has been shown in Section 4 that the existence of 
only small fractions f„ ~ 10~ 5 of protons in bound states, 
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FIG. 9: 4 He bound state fraction ft for a CHAMP BBN 
model (A) with Q. x h 2 = 0.1, t x = 3 x 10 4 s, and f E M = 1 
(the two curves on the left), and 1 H bound state fraction 
ft for a CHAMP BBN model (B) with fix h 2 = 5 x 10~ 3 , 
tx = 3 x 10 6 s, and /em = 1 (the two curves on the right). 
Solid (red) curves show ft when photodisintegration of bound 
states due to electromagnetic energy release during the X- 
decay is included, whereas dashed (green) curves show re- 
sults when this process is neglected. The resultant 6 Li yield 
in model (A) is ~ 10 times lower than when photodisin- 
tegration is excluded. Similarly, the resultant 6 Li yield in 
model (B) is ~ 100 times larger than without photodisin- 
tegration. CHAMP-exchange reactions have not been taken 
into account. 



TABLE IV: Rates for CHAMP-exchange reactions computed 
in the Born approximation. 



No. 




(AX) + B 


-> C + X 




rate [cm 3 s 


- 1 ] 


14 


OH 


-X-) + 2 H^ 


( 2 H-X") + 


: H 


8.8 x 10" 


15 


15 


Ch 


-A") + 3 H-> 


( 3 H-X") + 


4 H 


1.4 x 10" 


15 


16 


( 2 H 


-X~) + 3 H^ 


( 3 H-X") + 


2 H 


1.0 x 10" 


14 


17 


OH- 


X~) + 4 He^ 


(*He-X~) 4 


- *H 


3.6 x 10" 


17 


18 


( 2 H- 


X-) + 4 He^ 


( 4 He-A") 4 


- 2 H 


2.9 x 10" 


16 


19 


( 3 H- 


X~) + 4 He^ 


( 4 He-X") 4 


- 3 H 


8.0 x 10" 


16 



forming below T < 1 keV, may efficiently destroy again 
any priorly synthesized 6 Li and 7 Be. In Section 5 it has 
been seen that the efficiency of this destruction may be 
significanlty reduced when non-thermal photodestruction 
of bound states is taken into account. In this section, 
a further important process reducing late-time 6 Li and 
9 Be destruction is discussed. CHAMPs in bound states 
may exothermically transfer to heavier nuclei of equal 
or higher charge. In particular, ( 1 H— X~) bound states 
could be removed by the ( 1 H-A") + 4 Hc^ ( 4 He— A") 
+ 1 H charge exchange process. Charge exchange reac- 
tions turn out to be very important. In Table 4 the most 
important of these processes are presented. Rates for 
these processes were calculated in a very similiar way, i.e. 
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FIG. 10: As Fig. [5] but with charge exchange reactions in- 
cluded. 
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FIG. 11: As Fig. [7] but with charge exchange reactions in- 
cluded. 



in the Born approximation, to those of nuclear reactions 
involving bound states, as presented in Appendix A. Here 
the dipole (quadrupole) operators Eq. (|A4p (Eq. IA5|) 
replaced by the electromagnetic potential between the 
bound state and the heavier nucleus. The same argu- 
ments as presented in Section 3 apply concerning the fail- 
ure of the Born approximation. In particular, rates given 
in Table 4 should be only considered as benchmarks, with 
the true rates possibly deviating significantly. 

Fig. [TUl shows bound state fractions for the same model 
as that shown in Fig. 03 but now with CHAMP exchange 
reactions included (photodisintegration of bound states 
is neglected). From the comparison of these two figures 
it is evident that whereas bound state fractions of X H 
in the absence of exchange reactions reach levels close 
t° fp ~ 10~ 3 , they are two orders of magnitude below 
when exchange reactions are present. This is mostly due 
to the OH- A") + 4 He^ ( 4 He— A") + X H reaction. A 
for the final BBN yield almost equally important change 
is the elevated 2 H (and 3 H) bound state fraction when 
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the reactions in Table 4 are included. Though most 1 H 
exchange their CHAMPs with 4 He, due to the large 4 He 
abundance, a large fraction ~ 1 of 2 H enter bound states 
by capture of CHAMPs from protons as well. The 2 H 
bound state fraction in Fig. [10] (as well as Fig. [TTj) is 
only small /|, <C 1, simply because once a 2 H (and 3 H) 
bound state has formed, its life time against destruction 
by reactions shown in Table 3, is very short. In other 
words, essentially each 2 H which enters a bound state 
will be subsequently destroyed, leading to the produc- 
tion of 3 He and 4 He. This will have important conse- 
quences for bounds on CHAMPs at larger CHAMP den- 
sity, since either the lower bound on 2 H or the upper 
bound on 3 He/ 2 H may be violated. Fig. [11] shows the 
abundance evolution corresponding to Fig. 1101 and is the 
equivalent to Fig. [7] but now with exchange reactions 
switched on. Several trends are visible: With charge 
exchange reactions the 2 H/H ratio has fallen below the 
observational lower limit, i.e. 7.4 x 10" 6 compared to 
2 x 10~ 5 in Fig. [7] the final 7 Li/H ratio is larger, i.e. 
9.5 x 10" 10 compared to 2.7 x 10~ 10 , and the 6 Li/H ratio 
is much larger, i.e. 3.3 x 10~ 10 compared to ~ 4 x 10 -14 . 
Here 7 Li is larger due to reduced 7 Be( 4 H— X~, A"") 8 B 
and enhanced 4 He( 3 H— X~, A~) 7 Li efficiencies, and 6 Li 
is larger due to reduced 6 Li( 1 H— A", A") 4 He + 3 He and 
enhanced 4 He( 2 H— A", A") 6 Li reactions. 

When doing bound-state BBN computations with re- 
actions on Z = 1 bound states as well as CHAMP ex- 
change reactions included often very counter-intuitive 
results are obtained. As only one example, when the 
( 1 H— A") + 2 H rate is increased the 6 Li (and 7 Be) abun- 
dance may be reduced drastically. This is not what is ex- 
pected since a lower ( X H— A") and higher ( 2 H— X~) frac- 
tion ought to lead to a higher 6 Li abundance via enhanced 
4 He( 2 H-A-, A-) 6 Li and reduced 6 Li( 1 H-A-, A-) 4 He 
+ 3 He. Nevertheless, this is not what happens, due to a 
higher ( 2 H— A") fraction more 2 H is destroyed initially, 
rendering the 4 He( 2 H— A", A") 6 Li less effective at late 
times. Since the final abundance yield is given by the bal- 
ance of the still fast processes of 4 He( 2 H— X~, A") 6 Li 
production and 6 Li( 1 H-A", A~) 4 Hc + 3 He destruc- 
tion at late times less 6 Li results. Due to a lower 
6 Li( 1 H— X~, A~) 7 Be efficiency less 7 Be results. Late- 
time bound-state BBN is very non-linear requiring full 
numerical integration up to late times to obtain reliable 
predictions. 
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FIG. 12: Abundance yields of 2 H/H solid (red) 7 Li/H dashed 
(green), and 6 Li/H dotted (blue) as a function of CHAMP-to- 
baryon ratio Y x for a model with r x = 10 12 sec and excluding 
electromagnetic- and hadronic- energy injection. 



TABLE V: Three realizations of models which fulfill con- 
straints on light-element abundances and reconcile predicted 
with observed 7 Li/ 1 H and 6 Li/ 1 H ratios . Shown are the 
CHAMP-to-baryon ratio as well as a list of reaction numbers 
and the respective factors by which these reactions rates have 
been multiplied with respect to the (unreliable) estimates in 
the Born approximation. All models have r x = 10 12 s and 
electromagnetic or hadronic energy injection has not been 
taken into account, corresponding to an invisible or almost 
mass-degenerate decay. Abundance yields in these models 
are shown in Table IVT1 



Model 


Y x 


Reactions modified 


A 


4.3 x 10" 4 


#4: 0.1 #7: 2. #14: 0.1 #17: 0.3 


B 


4.3 x 10" 4 


#4: 0.3 #9: 0.3 #17: 0.1 #18: 30. 


C 


4.3 x 10~ 6 


#5: 3. #7: 30. #14: 0.03 #17: 0.03 



may be solved at once in the presence of a decaying par- 
ticle with decay time t x « 1000 s. This is possible in ei- 
ther case, a charged relic or a neutral relic. In subsequent 
sections it has been seen that late-time nucleosynthesis 
in the presence of charged weak-scale mass particles may 
lead to orders-of-magnitude modifications of the 6 Li, 7 Li, 
(and 2 H) abundances. It would be interesting to know 



VII. SOLUTIONS TO THE 6 LI AND 7 LI 
PROBLEMS DUE TO BOUND-STATE BBN FOR 
LONG-LIVED r^>10 6 SEC CHAMPS ? 



TABLE VI: The corresponding abundance yields resulting in 
the models shown in Table FVl 



In Section 2 priorly proposed solutions to the 7 Li 
overabundance and 6 Li underabundance resulting within 
BBN in the presence of (relatively) short-lived CHAMPs 
have been discussed. Notwithstanding possible astro- 
physical explanations of these deviations between theory 
and observation, it has been shown that both problems 



Model 


2 H/H 


7 Li/H 


6 Li/H 


A 


2.6 x 10" 5 


2.8 x 10" 


10 


9.3 x 10" 


12 


B 


2.4 x 10" 5 


2.3 x 10" 


10 


3.7 x 10" 


11 


C 


2.7 x 10" 5 


1.5 x 10" 


10 


3.2 x 10" 


11 



11 




FIG. 13: As Fig. [13] but for reaction rates as in Model C 
shown in Table IVT1 

if CHAMPs with long life times t x > 10 6 s may reconcile 
the 6 Li and 7 Li discrepancies. 

In Fig. rjJ2] abundance yields for t x = 10 12 s and vary- 
ing Y x axe shown. Here reaction rates in the Born 
approximation were adopted and electromagnetic- and 
hadronic- energy injection due to the X decay was ne- 
glected, corresponding to, for example, an invisible decay 
or a decay to a neutral daughter particle almost degen- 
erate in mass with the CHAMP. The model also approx- 
imates well the case of no decay, i.e. a stable CHAMP. 
At low CHAMP-to-baryon ratio Y x only 6 Li is modi- 
fied. Here most of the 6 Li is synthesized not at early 
times due to 2 H( 4 He— X~ , X~) 6 Li but rather at late 
times due to 4 He( 2 H-A~, A~) 6 Li. A small CHAMP 
density may therefore easily account for 6 Li in Pop II 
stars. When Y x increases to 10 -8 too much 6 Li is synthe- 
sized. For larger Y x > 10~ 3 6 Li destruction due to a high 
( 1 H-A' _ ) fraction reduces 6 Li again to observationally 
friendly levels. However, such models are then ruled out 
by 7 Li overproduction and 2 H underproduction, due to 
high ( 2 H-A~) and ( 3 H-A"~) fractions, with 7 Li produced 
by 3 H( 4 He— X~, A~) 7 Li and 2 H destroyed by reactions 
given in Table 3. When the decay is electromagnetic or 
hadronic, with a large fraction Jem ~ 1 of rest mass 
of X converted to electromagnctically interacting parti- 
cles such high Y x should in any case be ruled out due to 
elevated 3 He/ 2 H-ratios (cf. Ref. 

Nevertheless, significant uncertainties exist due to the 
uncertainties in the bound-state nuclear reactions and 
charge exchange reactions. In Tables [V] and I VII three 
(somewhat randomly chosen) models which do solve the 
6 Li and 7 Li problems are shown. Here a number of re- 
action rates were scaled up (or down) from the Born ap- 
proximation in order to arrive at an observationally sat- 
isfying result. It is seen that even at low Y x such models 
may exist, depending on the exact magnitude of rates 
for a variety of reactions. It is also seen that, when go- 
ing to lower Y x , rates have to deviate more drastically 
from the Born approximation in order to solve the 6 Li 



TABLE VII: Adopted values for ff ui for the different reac- 
tions varied in the Monte-Carlo analysis (see text for details). 
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FIG. 14: Probability in the CHAMP-to-baryon Y x - CHAMP 
life time r x parameter space at large r x that simultaneous 
solutions to the 7 Li and 6 Li (dark-blue) or only 7 Li (light- 
blue) problems exist. The points indicate 1 — 5% probability, 
whereas white areas had less than 1% of all randomly chosen 
reaction rates in the Monte-Carlo analysis result in 6 Li and 
7 Li (or 7 Li only) solutions. No electromagnetic- or hadronic- 
energy release has been taken into account. Areas above the 
red line would be ruled out due to electromagnetic cascade 
nucleosythesis under the assumption that /em—Sx 10~ 2 of 
the rest mass of the CHAMP is converted to electromagnti- 
cally interacting energy. See text for further detail. 

and Li problems. The corresponding abundance yields 
for Model C, where at low Y x observationally satisfying 
results are obtained, are shown in Fig. 1131 The figure 
clearly indicates that parameter space for a reduction of 
7 Li and production of some 6 Li exists. 

In the absence of reliable estimates for reaction rates 
it is difficult to assess quantitatively if significant param- 
eter space for simultaneous solutions for the 6 Li and 7 Li 
discrepancies for late decaying tx ^ 10 6 sec CHAMPs ex- 
ist. In particular all nuclear reactions shown in Table |TT] 
and Table IIII1 as well as the charge exchange reactions 
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shown in Table IIV1 i.e. a total number of nineteen re- 
actions. Though all rates have been determined numer- 
ically in the Born approximation in this paper, as the 
Born approximation is likely to fail badly, results become 
uncertain. In order to still arrive at a reliable result one 
is thus forced to perform a Monte-Carlo analysis, vary- 
ing all ill-determined reaction rates within conservative 
ranges. This has been done in the present paper. In par- 
ticular, the Born approximation values of the rates shown 
in Figs.[3]and dj as well as given in Table HTll and llVi have 
been taken as benchmarks. For each reaction a random 
generator determined a factor /j with which the bench- 
mark rate was multiplied. These factors where generated 
with a probability distribution flat in logarithmic space, 
and between values 1/ f[ ut < fi < f[ ut - For the reaction- 
rate dependent conservatively chosen ff ut the reader is 
referred to Table rVTil For each point in parameter space, 
i.e. for Y x and t x , this procedure was repeated a 1000 
times in order to arrive with one thousand different ran- 
domly chosen sets for the 19 ill-determined reaction rates. 
For each realization of reaction rates an indpendent BBN 
calculation was then performed and compared to the ob- 
servational constraints. 



The results of this Monte-Carlo analysis are shown in 
Fig. [14] Here dark (dark-blue) area indicates the proba- 
bility that between 1% - 5% (i.e. 10-50) of all indepen- 
dent 1000 BBN calculations with randomly varied rates 
respect the abundance limits on other light elements (as 
given in Ref. [13) while fulfilling 7 Li/ 1 H< 2.5 x 10" 10 and 
0.66 > 6 Li/ 7 Li> 0.03. Similarly, light (light-blue) areas 
indicate the same, but with now only the 7 Li discrepancy 
solved (i.e. 6 Li/ 7 Li< 0.03) is acceptable). It is noted that 
in none of the parameter space a probability > 5% for 
6 Li + 7 Li (or only 7 Li) solving areas is found, indicat- 
ing that the reaction rate combinations which may yield 
such solutions are rather rare. The liklihood for such sce- 
narios is even further diminished when electromagnetic- 
and/or hadronic- energy injection due to the decay of 
the particle is considered. In fact, when /em ~ 1 ah of 
the parameter space shown in Fig. [14] capable of solving 
the 6 Li+ 7 Li problems simultaneously (though at a < 5% 
liklihood), would be completely eliminated. Only when 
Jem is rather small, some area remains. This is shown 
by the (red) line for Jem = 3 x 10~ 2 corresponding, for 
example, to the decay of a stau f to a tau and gravitino, 
with the gravitino only 10% lighter than the stau. The 
area above the line is ruled out by overproduction of the 
3 He/ 2 H ratio due to 4 He photodisintegration. On the 
other hand, not shown in Fig. [TJ] are areas where only 
the 6 Li abundance as observed in Pop II stars may be 
produced. These exist plentiful, and at high probability, 
in particular at lower Y x < 10 -5 . It thus seems unliklcy 
that CHAMPs with t x > 10 6 sec may resolve the 7 Li prob- 
lem, though they could possibly constitute the source for 
the observed 6 Li at low metallicity. 



VIII. CONCLUSIONS 

In summary, I have presented results of a very de- 
tailed study of BBN in the presence of negatively charged 
massive particles (CHAMPs). Such particles have been 
shown to form bound states with nuclei towards the end 
of a conventional BBN epoch [H, U H [H, H3, H, HI 
and may alter BBN yields due to the catalysm of nuclear 
reactions [23[ . The present analysis attempts to take into 
account of all relevant effects for making relatively pre- 
cise predictions of catalytic light-element nucleosynthesis 
for nuclei with A < 7, but excluding the formation of 
molecules. It includes numerical evaluations in the Born 
approximation of all key nuclear cross sections, where 
one of the nuclei is in a bound state. Bound-state recom- 
bination and photodisintegration cross sections are also 
determined numerically. Furthermore, three very impor- 
tant and priorly not treated effects for the CHAMP BBN 
at late times r > 10 5 sec are included: (a) rapid nuclear 
reactions including charge Z = 1 nuclei in bound states, 
(b) the photodisintegration of bound states due to 7- 
and x- rays generated during the decay of the CHAMPs, 
and (c) CHAMP-exchange reactions from a bound state 
within a lighter nucleus to a bound state within a heavier 
nucleus. Light element abundances and bound state frac- 
tions are computed without approximations. The effects 
of hadronic and electromagnetic cascades due to CHAMP 
disintegration on light element abundances are properly 
taken into account. 

The present detailed study reveals that bound- 
state BBN proceeds very differently than initially fore- 
casted [23, [24J . At low temperatures T < 1 keV, a large 
number ~ 20 of Coulomb-barrier unsupressed nuclear 
reactions and charge exchange reactions become oper- 
ative and are capable, in most of the parameter space, to 
change 6 Li, 7 Li, and 2 H abundances by orders of mag- 
nitude. Unfortunately, reaction rates for these processes 
are not well approximated by the Born approximation, 
such that for CHAMP life times t x > 10 5 sec one has to 
resort to a Monte-Carlo analysis. 

The purpose of this study is to investigate the potential 
of CHAMP BBN to resolve the current 6 Li and 7 Li dis- 
crepancies between standard BBN and observations. It 
is shown, that a priorly proposed simultaneous solution 
of the 6 Li and 7 Li problems with a relic particle decaying 



1000 sec 21 



is not very dependent on the de- 
caying relic being charged 35j or not, unless its hadronic 
branching ratio is well below Bh "Z 10~ 4 . A solution with 
Bh <C 10 -4 has, however, the advandtage to not change 
much the 2 H/ 1 H ratio from its respective standard BBN 
value. Since 6 Li and 7 Li may be rapidly destroyed at late 
times one generically expects further simultaneous solu- 
tions of the 6 Li and 7 Li problems for t x > 10 6 sec. Nev- 
ertheless, even given the current reaction rate uncertain- 
ties, a Monte-Carlo shows that only a very small fraction 
^ 5% of reaction rate combinations may lead to such so- 
lutions. Since such possible solutions occur at relatively 
high CHAMP-to-baryon ratio 3 x 10" 5 < Y x < 10~ 2 they 
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are further constrained by the effects of electromagnetic 
energy injection and possible 3 He/ 2 H overproduction, re- 
quiring the decay to be invisible, or mother and daughter 
particle to be somewhat degenerate in mass ~ 10%. On 
the other hand, CHAMPs may well be the source of the 
observed 6 Li at low metallicity. 

I acknowledge helpful discussions with and M. As- 
plund, S. Bailly, O. Kartavtsev, K. Kohri, A. Korn, 
G. Moultaka, M. Pospelov, J. Rafelski, G. Starkman, 
V. Tatischeff, and T. Yanagida. 



APPENDIX A: THERMONUCLEAR REACTIONS 
IN THE PRESENCE OF BOUND STATES IN 
THE BORN APPROXIMATION 

Consider the three-body system of nuclei A, B, and 
CHAMP X~. Since for weak scale mass CHAMPs and 
light nuclei Mx 3> Ma , Mb , it is an excellent approxima- 
tion to assume X~ to be at rest at the origin, effectively 
acting as an external potential whichs absorbs momen- 
tum but not energy. The Hamiltonian of the system is 
then given by 

H = ^M A r 2 A + \m B v b + V c {\r A -t b \) + 

tt , \ z az 2 Z B e 2 

VNUc(rA,r B ) , (Al) 

ta r B 

where ta, represent the position vectors of nuclei A 
and B, ta, tb their magnitudes, and Zag, Zbg their 
respective charges. In Eq. (|Alj) the first two terms repre- 
sent kinetic energies, the second and third term, Coulomb 
and nuclear potentials between A and B, and the last 
terms, the Coulomb potentials between the (assumed 
singly charged) CHAMP X~ and the nuclei. This Hamil- 
tonian will be split into a dominant contribution Hq and 
a perturbative contribution H± <C Hq. In a rearrange- 
ment reaction of the type (A — X~) + B — ► C + X~ , 
where C is a nuclear bound state between A and B, the 
unperturbed and perturbed Hamiltonians for initial and 
final states are different, i.e. H\ ^ H$, H\ ^ H[. In 
particular, whereas in the initial state the perturbation 
is best chosen as the nuclear attraction beween A and B, 
i.e. H\ = Vnuc an d Hq = H — H[, in the final state 
it will be the differential Coulomb force of A~ on the 
nuclear bound state C = (A — B). When initial and final 
states are chosen as eigenstates to Hq and Hq, respec- 
tively, standard methods show that, in the Born approx- 
imation the transition amplitude may be computed by 
either (i\H[ /) or (/|i?gi). The initial and final states 
are chosen as 

\i) = \$ {A -x-)(rA)) \®Coui(r B )) (A2) 
I/) = \®(A-B)(P)) |*c«,i(«0) (A3) 

where s and p are the A — B center of mass and rel- 
ative coordinates, respectively. Coulomb wave func- 



tions \$coui) and the A — X~ bound state wave func- 
tion $m— x-) were determined numerically with realistic 
charge distributions. The nuclear wave function $>(a-b) 
was parametrised by <I> = 2-\/7 5 /3pexp(— jp) with 7 ad- 
justed such that in the absence of X~ the correct exper- 
imentally determined cross section results. The pertur- 
bation H( was chosen as the first non- vanishing element 
in the expansion of the last two terms of Eq. (|A1[) in 
terms of relative coordinate p. For dipole transitions this 
results into 

H{ = ~{Z A R A + Z B R B )e 2S -^ (A4) 
whereas for quadrupole transitons 

H[ = -(ZARl + Z B Rl)e^l^^-\^) (A5) 

where R A = Mb I (Ma + M B ) and R B = -M A /(M A + 
Mb)- Rates were evaluated by numerical integration of 
the matrix elements (i\H{ f) employing Fermi's Golden 
rule 

where V is a normalization volume, v relative velocity, S 
the Delta-function, and 

dN f = ^L^pldpcdQc (A7) 

a measure of the final phase space for nucleus C. For the 
evaluation of the matrix elements, six-dimensional inte- 
grals over the coordinates of two nuclei could be analyt- 
ically reduced to three-dimensional integrals which were 
numerically evaluated. Similiar to Ref. [3lj ] I have not 
considered internal spin of the nuclei, except for the obvi- 
ous total angular momentum degeneracy factors. Finally 
cross sections a(E) were converted to S-factors S(E) . 
They are related by 

a(E) = (S(E)/E) exp(-G(E)), (A8) 

where E is center-of mass (CM) energy and exp(G) with 

= 2n(Z A - l)Z B ac 

VCM 

is the Coulomb repulsion factor. In the above vcm is 
the relative velocity (vqm ~ «b for bound states), and 
a, c fine structure constant and speed of light, respec- 
tively. For assumptions concerning the angular momen- 
tum of the final A-B nucleus, the number of multipolcs 
included in the calculation, and the assumed S-factor in 
the absence of bound states the reader is referred to Ta- 
ble 2. The determined S-factors were subsequently in- 
tegrated over a thermal distribution to derive thermal 
nuclear rates in the presence of bound states. 
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